Dorathy  A.  Stewart 

Research  Directorate 
US  Army  Missile  Laboratory 


19  August  1981 


QD.S.ARMY  MQSSOILS 

fRatdmtonm  Arsenal,  A/abama  35898 


Approved  for  public  release;  distribution  unlimited. 


82  01  06012 


SMI  FOFW  1021,  1  JUL  79  PREVIOUS  EDITION  18  OBSOLETE 


DISPOSITION  INSTRUCTIONS 

DESTROY  THIS  REPORT  WHEN  IT  IS  NO  LONGER  NEEDED.  DO  NOT 
RETURN  IT  TO  THE  ORIGINATOR- 


DISCLAIMER 

THE  FINDINGS  IN  THIS  REPORT  ARE  NOT  TO  BE  CONSTRUED  AS  AN 
OFFICIAL  DEPARTMENT  OF  THE  ARMY  POSITION  UNLESS  80  DESIG¬ 
NATED  BY  OTHER  AUTHORIZED  DOCUMENTS. 


TRADE  NAMES 

U8E  OF  TRADE  NAMES  OR  MANUFACTURERS  IN  THIS  REPORT  DOES 
NOT  CONSTITUTE  AN  OFFICIAL  INDORSEMENT  OR  APPROVAL  OF 
THE  USE  OF  SUCH  COMMERCIAL  HARDWARE  UR  SOFTWARE. 


-  « ■: '  %V-'~  -  = 


■  i  sac-^if  ’  •  ■  - : 


J  [;!■ 


5*  <• 


SECURITY  CLASSIFICATION  O*  THU  PAGE  (Wi»fi  ?>»(«  Kntorod) _ 

_ REPORT  DOCUMENTATION  PAGE _ beforePcompLetwgN?Sqrm 

b  REPORT  NUMBER  2,  COVT  ACCESSION  NO,  S.  RECIPIENT’S' CRT ALOG  HUMBER 

Tt-RR-8i-A  f]  pi  '■/d  /A  '9  *v  /!  G> 


"A  ,101  b  c\o> 


4.  TITLE  (and  Subtttl*) 

A  Survey  of  Atmospheric  Turbulence 
Characteristics ^ 

7 7  AUTHORf*;  “ 


5.  TYPE  OF  REPORT  ft  PERIOD  COVERED 


Technical  Report 


8.  PERFORMING  ORG.  REPORT  NUMBER 


e.  CONTRACT  OR  GRANT  NUMBERS 


Dorathy  A.  Stewart 


»■  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Commander,  US  Army  Missile  Command 

ATTN:  DRSMI-RRA 

Redstone  Arsenal,  AL  35898 

II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Commander,  US  Army  Missile  Command 

ATTN:  DRSMI-RPT 

Ked stone  Arsenal,  AL  35898 


10.  PROGRAM  ELEMENT,  PROJECT,  TASK 
AREA  Ik  WORK  UNIT  NUMBERS 


12.  REPORT  DATE 

19  August  1981 


(JISOIIJ.  IU  1  13  NUMBER  OF  PAGES 

Redstone  Arsenal,  AL  35898  33 

MONITORING  AGENCY  NAME  ft  AODRESSfi/  dltloronl  from  Controlling  OWco)  15.  SECURITY  CLASS,  (nl  thl.  report) 


Unclassified 

15*.  DECLASSIFICATION/ DOWNGRADING 

schedule 


16.  DISTRIBUTION  STATEMENT  (of  thl.  R.port) 

Approved  for  public  release;  distribution  unlimited. 


17.  DISTRIBUTION  STATEMENT  (of  thm  .b.tr.ct  anfarad  In  Block  30,  II  dlllortnt  Irom  R.port) 


...o 

«?'  4  v  % 


1$.  SUPPLEMENTARY  NOTES 


Y  's*  at 


19.  KEY  WORDS  (Conf/nu#  on  rmvno  aidm  If  n«c wnd  identity  by  block  number; 

Turbulence 

Wind 

Planetary  Boundary  Layer 


1 20.  ABSTfnAqT  fCot»</mi»  on*  rmrM  eftf v  f/  necMierf  mad  Mmntlty  by  block  numbmr) 


This  report  reviews  the  literature  on  atmospheric  turbulence  in  the 
planetary  bounda.  ?  layer,  with  the  main  emphasis  on  the  lowest  tens  of  meters. 
A  table  summarizes  26  references  which  discuss  variations  in  intensity  of 
turbulence.  This  intensity  tends  to  be  greater  over  rough  surfaces  than  over 
smooth  surfaces.  Intensity  of  turbulence  normally  decreases  when  either 
altitude  or  atmospheric  stability  increases.  High  wind  speeds  are  more  likely 
to  be  associated  with  nearly  neutral  stability  than  with  very  stable  or  very 


DD  1  jaJTVs  1473  EDITION  OF  »  NOV  88  IS  OBSOLETE 


_ UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  (Wran  Dolt  E.iforerf) 

L/J  /??>%■ 


unstable  conditions.  This  report  also  includes  a  discussion  of  the 
autocorrelation  function  for  the  three  components  of  the  wind  vector. 


_  UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAOEft»Ti»n  Cat*  Bntirtd) 


rt.-i'Sfi  •bsnpan^'qtfrr 


ACKNOWLEDGMENTS 


The  author  expresses  appreciation  to  Dr.  Oskar  M.  Essenwanger 
for  encouragement  and  helpful  suggestions  during  this  project;  to 
Mr.  Richard  E.  Dickson  and  Dr.  George  Fichtl  for  suggesting  a  few 

references;  and  to  Mrs.  Louise  H.  Cooksey  for  her  clerical  assist- 

/ 

ance.  /  / 


Accession 

OTIS 
DT1C  T«3 

Unonno,.v>---'“d 


Br . .  — - - 

Pl.strib'ai'vn/ 

‘"Avail- /.-my  coaoa_ 

-  •  .  .  i  nr'Vor 

Ti  -4  *■  ^  i  .  0  V C  ^ 


TABLE  OP  CONTENTS 


Page 


I.  Introduction . . . . .  3 

II.  Intensity  of  Turbulence  ...................  4 

III.  Autocorrelations  .........  .  ......  15 

IV.  The  Ekman  Layer . . . .  18 

V.  Summary  . . 20 


LIST  OF  TABLES 


Table  Title  Page 

I.  Intensities  of  Turbulence  at  Test  Area  5  on 

Redstone  Arsenal  (u  *  m/seo)  .  5 

II.  Summary  of  Longitudinal  Intensities  of  Turbulence  .  8 


..*• HI' 


I. 


INTRODUCTION 


j 


-•  !»««es%M»*«JWaSM3 


i 


In  studies  of  atmospheric  turbulence,  the  positive  x  direction  ia  | 

usually  defined  as  the  direction  of  the  mean  horizontal  wind,  and  the  com¬ 
ponent  of  the  wind  in  the  x  direction  ia  the  u  component.  The  time  average  of 
u  is  denoted  by  u,  and  u'  *  u  -  u.  The  other  horizontal  coordinate  ia  y,  and 
the  component  of  the  wind  in  the  y  direction  is  v.  By  definition  the  time 
average  of  v,  v,  ia  zero,  and  v*  *  v.  The  vertical  coordinate  ia  z,  and  the 
xyz  system  ia  a  right-handed  rectangular  coordinate  system.  The  vertical  com¬ 
ponent  of  the  wind  ia  w.  The  quantities  u*  and  v*  are  often  referred  to  as 
the  longitudinal  and  the  lateral  fluctuations  of  velocity,  respectively.  The 

corresponding  standazd  deviations  are  au  and  ov,  where  ou  «*  [ ( u '  j V?  and  av  * 

[(7^]^.  The  longitudinal  intensity  of  turbulence  is  defined  aa  iu  *  ou/u, 

and  the  lateral  intensity  is  iv  *=  cv/u. 

In  approximately  the  lowest  50m  the  horizontal  wind  stresses  are 
nearly  constant  with  height,  and  the  wind  direction  ia  also  nearly  constant 
[Hess,  195u].  Then  the  following  equation  can  be  derived: 


5iT  _  u*  ( 1  ) 

9z  kz  * 


where  u#  ia  called  the  friction  velocity  and  k  is  the  von  Karman  constant. 
[Busch,  197?]  discuases  some  different  estimates  of  the  von  Karman  constant, 
but  it  is  now  most  commonly  taken  to  be  0.4*  [Panofsky,  197'?]  shows  that 
the  friction  velocity  at  100m  is  typically  only  about  10  per  cent  less  than 
the  surface  value  in  mid-latitudes. 

Integration  of  Equation  (1)  yields  u  =  (u#/k)Jlnz  plus  a  constant, 
but  the  logarithm  of  zero  is  undefined.  Therefore,  a  roughness  parameter,  z0, 
is  introduced  to  obtain 

*  "  (to  *o  ’  *)  '  (2) 

where  ^  is  a  stability  parameter.  This  wind  law  does  not  apply  below  z  *  z0. 

Tables  of  roughness  parameter  appear  in  [Hess,  1959]  and  [Frost  et 
al.,  1978].  Over  ice  z0  nay  be  less  than  0.01cm.  In  high  grass  or  wheat 
s0  is  typioally  a  few  centimeters*  In  a  forest  the  roughness  parameter  is  a 
fraotion  of  a  meter  to  one  meter.  In  a  city  the  roughness  parameter  is  nor¬ 
mally  1-4m. 

'?  is  zero  under  conditions  of  neutral  stability  and  the  simple 

i  .  logarithmic  wind  law,  u  -  (u*/k)  Jin  (z/z0)  is  obtained.  Neutral  stability 

\  exists  when  the  lapse  rate  of  temperature  is  equal  to  the  adiabatic  lapse 

rate.  Under  this  condition  a  parcel  of  air  displaced  vertically  experiences 
no  buoyant  acceleration. 
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When  the  temperature  deoreaaee  with  altitude  move  rapidly  than  the 
adiabatio  lapse  rate,  the  atmosphere  is  unstable  and  4>  is  positive . 

[Blackadar  et  al,  1974]  contains  information  for  estimating  41,  which  may  be 
greater  than  1.7  for  a  very  unstable  atmosphere. 

In  a  stable  atmosphere  4/  is  negative  and  can  be  expected  to  have  a 
magnitude  less  than  1.0. 

[Pasquill,  196l]  devised  a  classification  scheme  for  estimating 
atmospheric  stability  from  surfaoe  (10m)  wind  speeds  and  amount,  of  heating  or 
cooling.  Pasquill  stability  olasses  range  from  A  through  P,  where  A  is  very 
unstable,  D  is  neutral,  and  P  is  the  most  stable  class.  According  to 
Pasquill' s  classification,  as  wind  speeds  increase  the  atmosphere  approaches 
neutral  stability.  When  speeds  are  greater  than  6m/sec,  conditions  are 
neutral  unless  there  is  strong  insolation,  in  which  case  the  atmosphere  may  be 
slightly  unstable.  At  night,  cloud  cover  less  than  50  percent  permits  enough 
cooling  of  the  surface  that  the  atmosphere  is  in  the  most  stable  P  Pasquill 
category  for  wind  speeds  of  3m/sec  or  less.  Whenever  there  is  a  heavy  over¬ 
cast  of  clouds,  the  atmosphere  is  neutral  regardless  of  wind  speed  or  time  of 
day  or  night.  During  the  day  the  amount  of  insolation  depends  upon  sun  angle 
and  cloud  amount  and  type.  [Luna  and  Churoh,  1972]  outline  details  of  a  pro¬ 
cedure  for  determining  insolation  from  standard  meteorological  observations. 

This  report  contains  an  extensive  discussion  of  variations  of  inten¬ 
sity  of  turbulence  with  atmospheric  conditions,  surface  roughness,  and  alti¬ 
tude.  There  is  also  a  discussion  of  the  autocorrelation  function.  Finally, 
other  conditions  affecting  turbulence  are  discussed. 

II.  INTENSITY  OF  TURBULENCE 

This  section  presents  a  discussion  of  some  measurements  of  intensity 
of  turbulence  made  at  Redstone  Arsenal,  Alabama,  and  a  survey  of  the  litera¬ 
ture  describing  results  obtained  by  others. 

Table  I  contains  u  and  the  longitudinal  and  lateral  intensities  of 
turbulence  measured  in  August  1973  under  unstable  atmospheric  conditions.  The 
site  was  the  Army  Gas  Dynamics  Laser  Range.  It  consisted  of  a  grass-covered 
plot,  61m  wide  and  655m  long,  with  trees  lining  each  side.  Five  towers  were 
located  137m  apart,  and  a  sixth  was  9* 1m  from  the  middle  tower  on  a  line  per¬ 
pendicular  to  the  line  of  the  other  five  towers.  More  information  can  be 
found  in  [Stewart,  1975]* 

Both  longitudinal  and  lateral  intensities  of  turbulence  vary  with 
height.  At  10m,  iu  varies  from  0.25  to  0.98  and  has  a  mean  of  0*56,  while 
iv  has  a  muoh  larger  variation  from  0.16  to  1 . 4 1  and  a  mean  of  0.53*  If  a  mean 
is  taken  of  the  individual  ratios  iv/iu,  0*94  is  obtained  at  10m.  At  6m 
iv/iu  *>as  a  mean  of  1.04.  The  mean  iv  of  0.61  is  also  greater  than  the  mean 

iu,  which  is  0.59*  At  2m  the  mean  of  the  ratios  iv/iu  is  1*10*  At  2m 

iv  varies  from  0.21  to  3*33  and  has  a  mean  of  0.85.  At  2m  iu  also  has  a  wide 

variation  from  0.27  to  2.09,  and  the  mean  is  0.77.  Note  that  iv/iu  “  0y/0u. 

It  follows  that  under  the  conditions  of  these  measurements  the  standard 
deviation  of  the  lateral  component  of  the  wind  is  the  same  order  of  magnitude 
as  the  standard  deviation  of  the  longitudinal  component. 
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TABLE  I.  INTENSITIES  OF  TURBULENCE  AT  TEST  AREA  5  ON  REDSTONE  ARSENAL  (tr  -  a/s ec) 


TABLE  I.  (Concluded) 
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Thera  is  no  agreement  among  diverse  seta  of  observations  on  the 
magnitude  of  ov/°u*  This  probably  represents  real  changes  in  atmospheric  con¬ 
ditions*  For  example,  [Swanson  and  Cramer's,  196$]  Table  6  contains  mean 
values  of  ov/°u  at  and  th®  magnitude  of  this  ratio  varies  from  0.44  to 

1.03.  [Mayer's,  1981]  Table  3  summarises  4  sets  of  measurements  made  in  a 
spruoe  forest,  and  here  the  :«tio  av/ou  varies  from  0.6$  to  0.86. 

There  is  evidence  that  the  ratio  ov/ou  depends  upon  stability.  For 
example,  [Wyngaard  and  Clifford,  1977]  inolude  in  their  Table  1  mean  values  of 
(oy/tr)2  and  (<jy/tf)2  for  5.7m  over  a  flat  uniform  Kansas  plain.  For  very 
unstable  ooaditions  the  mean  ov  divided  by  the  mean  cu  is  1*32.  For  moder¬ 
ately  unstable  conditions  this  ratio  is  0.95,  but  when  the  atmosphere  is  stable 
the  ratio  is  0.74.  From  Table  2  of  [Panofsky  et  al.,  1978]  one  obtains 
magnitudes  of  the  ratio  cy/cu  of  0.72  to  1.18  for  unstable  conditions  at  a 
height  of  2m.  [Champagne  et  al.,  1977]  give  oY  and  cu  for  4m  above  flat  fern 
land  during  unstable  conditions.  For  all  of  the  4  sets  of  data  in  their  Table  1 
the  ratio  ov/ou  is  greater  than  1.0  and  in  one  test  run  cv  was  more  than 
50  percent  greater  then  cu.  On  the  other  hand,  [Ariel  and  Nadeshina,  1976]  sum¬ 
marised  field  measurements  for  neutral  stratification,  and  ov/cu  varied  from 
0.74  to  0.89.  [skibin,  1972]  obtained  cv/o u  of  0.48  in  an  experiment  where 
stability  was  considered  neutral  according  to  its  Pasquill  category,  and 
Cy/oy  was  0.26  for  slightly  unstable  conditions.  For  more  unstable  conditions 
Skibin  obtained  cv/cu  from  0.27  to  0.37.  The  problem  of  making  these  data 
consistent  with  other  investigations  may  be  the  inexactness  of  the  Pasquill 
method  [Luna  and  Church,  1972]. 

The  relative  magnitudes  of  cv  and  cu  uIbo  depend  upon  height. 

[Frost  et  al.,  1978]  consider  a  neutrally  stable  atmosphere,  and  from  their 
equation  4.26,  which  applies  at  a  height  of  10m,  oY/cru  -  0.64  is  obtained. 

They  then  proceed  to  give  a  typioal  example  where  this  ratio  increases  to  1.00 
at  600m,  above  which  cv  =  au.  According  to  [Dickson  and  Angell's,  1968] 

Figure  5  ov  ■  ou  at  2km,  but  av  is  from  0.7  ou  to  0.9  cu  0*5km.  [Bowne  and 
Ball,  1970]  list  means  of  Oy/u*  and  Cu/u*  taken  from  different  stabilities  for 
two  levels  in  a  rural  location  and  two  levels  in  an  urban  location.  Outside 
the  city  the  mean  cv  divided  by  the  mean  ou  is  0.76  at  12.2m  and  0.97  at  61m. 
Inside  the  city  the  ratio  is  0.76  at  15. 3m  and  0.60  at  53.3m.  [Bradley's,  1980] 
Table  3  describes  measurements  at  the  crest  of  a  170-m  hill;  these  data  do  not 
show  a  consistent  change  with  height.  cv/cu  ranges  from  0.65  to  0.74  at  9n> 
and  from  0*77  to  0.82  at  16m.  At  25m  the  ratios  vary  from  1.09  to  1.17; 
however,  at  87m  the  magnitudes  of  ov/ou  are  lower  than  at  25m  and  range  from 
0.71  to  1.02. 

Table  II  summarizes  some  of  the  literature  describing  measurements 
of  iu  and  ou/u#.  The  magnitude  of  iu  varios  from  less  than  0.1  to  more  than 
1.0,  but  the  most  typical  values  are  between  0.1  and  0.4*  The  values  of 
a^/ u*  are  commonly  between  1.5  and  3* 5* 

The  intensity  of  turbulence  normally  decreases  as  stability 
increases.  This  is  illustrated  in  the  work  of  [Swanson  and  Cramer,  1965]  who 
made  measurements  at  White  Sands  Missile  Range.  Their  tower  stood  on  a  smooth 
plot  surrounded  by  ground  containing  small,  uniformly  distributed  sand  dunes. 
They  measured  temperature,  wind  speed,  and  wind  direction  at  nine  levels  from 
4.6  to  62.0m.  Only  observation  periods  with  regular  temperature  profiles  and 
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TABLE  II.  SUMMARY  OF  LONGITUDINAL  INTENSITIES  OF  TURBULENCE 
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"ABLE  II.  (Cantlnued) 


4.6-m  wind  speeds  of  at  least  4  mi/hr( 1 .8m/sec)  were  included  in  the  data 
base.  They  used  the  rather  unusual  definition  that  conditions  were  neutral  if 
the  magnitude  of  the  temperature  difference  between  4*6m  and  62.0m  was  not 
greater  than  0.56°C.  On  the  basis  of  824  observation  periods  Swanson  and 
Cramer  found  that  when  the  atmosphere  was  less  stable  than  neutral,  the  inten¬ 
sity  of  turbulence  was  greater  than  when  the  atmosphere  was  stable. 

Intensities  were  intermediate  during  neutral  stability  at  each  of  the  nine 
levels.  For  example,  at  62.0m  for  wind  speeds  from  1.8  to  5*1  m/sec,  iu  * 

0.30  for  unstable  conditions,  iu  =0.16  for  neutral  conditions,  and  iu  ■  0.13 
for  stable  conditions.  At  4«6m  the  corresponding  intensities  are  0.42,  0.27, 
and  0.25,  respectively.  At  higher  wind  speeds  turbulence  intensities  are 
smaller.  For  wind  speeds  greater  than  4«46m/sec  under  unstable  conditions  the 
longitudinal  intensities  of  turbulence  are  0.22  and  0.17  at  4*6  and  62.0m, 
respectively.  Under  stable  conditions  the  corresponding  magnitudes  of  iu  are 
0.19  and  0.10. 

[Skibin,  1972']  includes  a  table  of  five  experimental  sets  of  obser¬ 
vations  made  in  connection  with  an  atmospheric  dispersion  study.  One  of  these 
was  under  Pasquill  stability  category  D,  or  neutral,  and  had  a  value  of  0.356 
for  cru/u.  In  a  slightly  unstable  case  cu/u  was  0.414*  For  more  unstable 
cases,  longitudinal  turbulence  intensities  were  0.309,  0.403*  and  0.852. 

[Wyngaard  and  Clifford,  1977]  summarize  earlier  work  in  their  Table  1 
in  which  they  list  the  mean  values  of  cu^/u^  over  a  flat,  uniform  Kansas 
plain.  Their  data  show  that  for  very  unstable  conditions  ou/u  is  0.22  and  is 
0.20  for  moderately  unstable  conditions.  In  a  moderately  stable  atmosphere 
ou/u  is  only  0.16. 

...  1 

In  [Binkowski* s,  1 97 ]  Figure  4  many  cases  of  ou/u#  are  plotted 
graphically  as  a  function  of  stability.  Near  neutral  conditions  the  values  of 
ou/ are  mostly  near  2.5,  and  under  stable  conditions  the  mean  is  nearer  2.2. 
For  an  unstable  atmosphere  the  mean  ou/u*  becomes  larger  as  the  atmosphere 
becomes  less  stable,  and  is  near  4.0  for  a  very  unstable  atmosphere. 

[Grimm,  1 97 1 ]  examined  15  cases  which  were  either  neutral  or 
slightly  stable,  where  z0  =  1.1cm,  and  obtained  a  mean  ou/u#  of  2.36  for 
heights  from  8  to  32m.  In  44  unstable  cases,  where  z0  =  2.6cm,  the  mean  Ou/u* 
was  3*31  * 

The  intensity  of  turbulence  depends  upon  the  roughness  of  the 
underlying  surface,  as  well  as  upon  the  stability.  Usually  the  intensity 
increases  as  stability  decreases  or  as  the  roughness  of  the  underlying  surface 
increases,  [llanna,  1981  ]  points  out  that  an  exception  to  this  rule  may  occur 
in  very  stable,  light  wind  conditions. 

[Bowne  and  Ball,  1970]  obtained  observations  of  turbulent  wind  fluc¬ 
tuations  on  a  tower  in  downtown  Fort  Wayne,  Indiana,  and  on  another  tower  in  a 
nearby  rural  setting.  The  roughness  was  larger  in  the  city,  and  the  urban 
heat  island  reduced  atmospheric  stability,  especially  at  lower  levels. 
Turbulence  was  more  intense  in  the  rougher  and  less  stable  urban  environment. 
The  lower  levels  on  the  urban  and  rural  towers  were  15*3  and  12.2m, 
respectively.  In  16  of  19  testa  au/il  at  the  lower  level  was  greater  on  the 
urban  than  on  the  rural  tower,  and  in  many  of  these  16  tests  the  difference 
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was  quite  large-.  The  upper  levels  on  the  urban  and  rural  towers  were  55*3  and 
61.0m,  respectively.  The  longitudinal  intensity  of  turbulence  in  the  urban 
location  was  greater  in  14  of  18  tests.  On  the  other  hand  when  one  considers 
°u/u*  one  finds  that  at  the  upper  level  the  mean  of  2.48  in  the  urban  location 
is  only  slightly  larger  than  the  2.42  for  the  rural  setting.  At  the  lower 
level  in  the  city  the  mean  cru/u#  of  4 •  1 6  was  much  larger  than  the  2.47 

measured  outside  the  city.  •'  .  .  i.  . 

[SethuRaman  and  Raynor,  1 980 3  compared  longitudinal  intensities  of 
turbulence  at  a  height  of  8m  over  the  Atlantic  Ocean,  151cm  from  Long  Island, 

New  York,  with  simultaneous  measurements  at  8m  above  the  beach.  As  with  pre¬ 
viously  discussed  studies  over  land,  intensit-.  of  turbulence  over  the  ocean 
decreases  as  stability  increases.  An  overall  average  cu/u  over  the  ocean  is 
near  0.09,  with  a  variation  from  approximately  0.02  to  0.19*  The  behavior  of 
the  ratio  of  longitudinal  intensity  of  turbulence  over  the  ocean  to  that  over 
the  beach  depends  upon  whether  the  flow  is  basically  onshore,  offshore,  or 
along  the  shore.  When  flow  is  offshore  the  intensity  of  turbulence  over  the 
ocean  is  approximately  the  same  as  that  over  land.  When  flow  is  along  the 
shore,  turbulence  intensity  over  the  ocean  is  about  half  that  over  land  when 
the  land  wind  speed  is  greater  than  6m/sec;  but  for  wind  speeds  less  than 
3m/sec  turbulence  intensity  over  the  ocean  is  greater  than  that  over  land. 

For  onshore  winds,  the  intensity  of  turbulence  over  land  and  water  is  about 
the  same  when  the  wind  speed  over  land  is  greater  than  12m/sec,  and  there  is  a 
minimum  ratio  of  oceanic- to- land  turbulence  of  0.5  near  10m/sec.  Below 
lOm/see  the  ratio  of  intensity  over  the  ocean  to  that  ever  land  increases  to 
almost  2  as  wind  speed  decreases  for  onshore  flow.  :  ! 

[tllen's  1968]  Table  4  summarises  measurements  made  near  Itchaca, 

New  York,  in  a  plantation  of  Japanese  larch  which  had  a  mean  height  of  10.40m. 
At  this  height  intensity  of  turbulence  was  0.47,  and  it  increased  to  0.57  at 
7*25m.  Variation  of  intensity  was  irregular  down  to  1.15m  where  it  was  0.51* 

[?ionco,  1972]  discusses  canopies  3uch  as  rice  paddies,  wheat 
fields,  and  forests.  Just  above  the  different  types  of  canopies  the  intensity 
of  turbulence  varies  from  0.28  to  0.47,  but  within  canopies  intensities  range 
from  0.32  to  0*84. 

['layer,  1981  ]  found  that  the  longitudinal  intensity  of  turbulence  at 
the  top  of  a  spruce  forest  varied  from  0.59  to  1.J0. 

In  the  previous  discussions  of  variation  of  intensity  of  turbulence 
with  surface  roughness  and  with  atmospheric  stability,  the  reader  may  have 
noticed  that  there  also  appeared  to  be  variations  with  altitude.  In  general, 

nu/7  can  be  expected  to  decrease  as  altitude  increases.  There  is  also  evi¬ 
dence  that  ou  decreases  with  altitude. 

[<ichtl  and  McVehil,  1970]  considered  a  large  amount  of  data  to 
develop  their  Table  2  in  which  the  ratio  aj  is  equal  to  2.227  under 

neutral  conditions  and  1.897  under  unstable  conditions,  where  u#Q  is  the  sur¬ 
face  friction  velocity.  Fiehtl  and  McVehil 's  Table  1  gives  5U  as  a  function 
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of  z  in  meters.  For  neutral  conditions,  By  ■  (z/18)“®*^,  and  for  unstable 
conditions,  By  -  (z/18)“°* *4.  Thus,  one  obtains  Oy/u^o  ■  2.227 
( a/1 Q) 31 5  for  neutral  conditions  and  cu/u#o  ■  1 .897(b/18)"®*^  under 
unstable  conditions. 

[Swanson  and  Cramer,  1965]  analyzed  both  the  longitudinal  and 
lateral  intensities  of  turbulence  on  a  62-m  meteorological  tower  over  a  two- 
year  period  at  White  Sands  Missile  Range.  Both  intensities  decreased  with 
height  in  all  thermal  stratifications.  Swanson  and  Cramer  found  that  the 
decrease  could  be  expressed  as  z  to  a  power  which  varied  from  -0.1  to  -0.3. 

The  magnitude  of  the  exponent  is  larger  for  more  stable  conditions.  They  also 
found  that  turbulent  intensities  at  all  heights  and  in  all  thermal  stratifies- 
;  tions  tended  to  be  inversely  proportional  to  the  mean  wind  speed. 

[DeLarrinaga,  1972]  tested  the  power  law  proposed  by  [Swanson  and 
Cramer,  1965]  on  wind  measurements  from  two  urban  sites  in  Liverpool.  A  cap- 
\  tive  balloon  was  used,  and  the  upper  anemometer  was  at  305m.  The  height  of  a 
lower  anemometer  was  varied.  These  urban  measurements  verified  that  ou/u 
decreased  as  z  increased.  DeLarrinaga  fitted  the  data  to  the  power  law 
described  by  Swanson  and  Cramer  and  obtained  exponents  from  -0.14  to  -0.36. 

[Bowne  and  Ball,  1970]  made  simultaneous  measurements  at  a  rural 
and  urban  site.  At  the  rural  site  wind  measurements  were  at  12.2m  and  61  «0m. 
At  the  urban  site  instruments  were  at  15»3m  and  53»3m.  At  the  urban  locations 
there  were  18  sets  of  data  where  both  upper  and  lower  values  of  <JU/{T  were 
available;  in  all  18  cases  the  intensity  of  the  lower  level  turbulence  wsb 
greater  than  the  intensity  of  the  upper  level  turbulence.  In  14  of  17  sets  of 
measurements  from  the  rural  site  the  intensity  of  turbulence  was  greater  at 
the  lower  level  than  at  the  upper  level. 


[Petit  et  al.,  1976]  show  a  plot  of  ou/u  within  and  above  a  forest 
in  their  Figure  5*  Tnere  is  some  irregularity  within  the  forest  to  approxi¬ 
mately  2m  above  tree  top  and  then  a  decrease  with  altitude. 

[Bradley,  1980]  made  wind  measurements  on  a  tower  placed  on  top  of  a 
170-m  hill  during  atmospheric  conditions  associated  with  neutral  stability. 

In  the  three  sample  cases  in  Bradley's  Figure  3,  ou^S  at  9m  is  approximtely  3 
times  Oy/iT  at  87m.  In  the  text  they  ,  define  an  intensity  of  turbulence  e.3 

[d/3)  (7*  +  T7?  +  ^)]ls  /u  and  consider  means  of /this  intensity  for  all  data 

i  ! 

0.326  for  9m;  0.168  for  25m;  and  0.120  for  87m. 

[Dickson  and  Angell,  1968]  contain  figures  of  <ju  and  cu/u  as  func¬ 
tions  of  height  to  2km,  and  both  quantities  fiecrease  with  increasing  altitude 
during  the  one  summer  of  data  plotted;  however,  du  only  decreases  slightly, 
and  ov  increases  slightly  with  increasing  altitude.  On  the  other  hand,  cu/u 
decreases  from  0.43  at  0.5km  to  0.22  at  2km  and  cv/ 3  decreases  from  0»33  to 
0.21. 
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[Duchene-Marullaz,  1975]  found  that  longitudinal  intensities  of  tur¬ 
bulence  in  a  suburban  area  near  Nantes,  France,  decreased  from  0«30  at  10m  to 
0.20  at  60m  for  SSW  and  SW  winds.  The  decrease  was  from  0.28  at  10m  to  0.17 
at  60m  for  westerly  winds. 

( 

III.  AUTOCORRELATIONS 

In  this  section,  a  detailed  discusnion  of  autocorrelation  functions 
is  prefaced  by  an  explanation  of  their  importance.  Equations  are  derived  for 
the  variance  of  the  difference  between  u(t+x)  and  u(t)  and  for  the  variance  of 
a  random  component  of  turbulence  which  is  uncorrelated  with  the  turbulent 
fluctuation  either  at  that  time  or  at  another  time.  Each  of  these  equations 
contains  an  autocorrelation  function  R(x)  which  is  the  correlation  between  the 
values  of  u  at  two  times  separated  by  the  time  interval  x. 


Let  Au  be  defined  by  t.he  equation 

Au  «  u  (t+x)  -  u(t) 


(3a) 


where  t  is  time.  If  "u  is  assumed  to  be  constant  throughout  the  time  period 
being  studied,  Equation  (3a)  can  he  rewritten  as 

\ 

Au  -  u’(t+x)  -  u’(t).  (3b) 

Both  sides  of  Equation  (3b)  may  be  squared  and  averaged  to  obtain 


(Au)2  -  [u’(t+x)];-  +(u’(t)]2  -  2  u’  (t+x)u’  (t) .  (4a) 

Each  of  the  first  two  terms  on  the  right  hand  side  of  Equation  (4a)  is  equal 
to  Ou2,  and  the  third  term  equals  2ou2(R(x)).  Thus  is  obtained 

(Au)2  -  2  ou2(l-R(x)).  (4b) 

When  i-O,  R(x)  -  1,  and  (Au)2  is  zero  as  expected.  When  x  is  very  large, 
the  autocorrelation  is  zero,  and  the  variance  of  Au  is  twice  the  varianoe  of 


It  is  sometimes  convenient  to  assume  a  linear  relationship  between 
u’(t+x)  and  u*(t)  and  to  write 

u* ( t+x)  -  u‘ (t)R(x)  +  u"(t)  (5a) 

where  u"  is  independent  of  u' .  [Hanna,  1979]  recommends  such  a  relationship 
when  the  coordinate  system  is  following  an  air  parcel,  but  it  can  also  be  used 
to  describe  behavior  at  a  point  which  is  stationary  relative  to  the  earth. 

One  can  rewrite  Equation  (5a)  as 

u"(t)  -  u'(t+x)-u'(t)R(x).  (5b) 
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If  squares  and  averages  are  r,nde  of  both  sides  of  Equation  (5b)  the  following 
is  obtained: 


ou"2  *  ou,2R2(x)  +  ou'2-2R(x)u'  (t)u'  (t+x) .  (6a) 


Since  u' (t)u' (t+x)  is  equal  to  R(x)ou'2,  Equation  (6a)  can  be  rewritten  as 


°u"2  *  au,2(1-R2(x)) .  (6b) 

It  is  obvious  that  for  very  large  values  of  x  the  variance  of  uM  is  equal  to 
-  the  variance  of  u' • 

'  Information  can  be  obtained  about  spatial  correlations  by  applying 

Taylor's  hypothesis  which  relates  spatial  correlations  to  time  lag  correla¬ 
tions.  This  frosen  field  or  frosen  turbulence  approximation  makes  the  substi- 

_  tution  x  ■  Ax/T i,  where  Ax  is  the  spatial  lag  in  the  direction  of  3.  It  is 

^assumed  that  the  turbulence  is  homogeneous  in  the  x  direction  and  stationary 
in  time  [Lumley  and  Panofsky,  19641  or  (Webster  and  Burling,  198j)* 

In  order  to  test  the  relationship  between  space  and  time  correlation 
;  functions,  [Cramer,  1959]  carefully  selected  six  experiments  from  Project 

Prairie  Grass.  In  each  one,  the  observed  wind  uirection  was  within  25  degrees 
of  the  longitudinal  axis  of  the  instrument  array.  In  Cramer's  Figure  11  for 
daytime  and  Figure  12  for  nighttime  experiments,  both  spatial  and  temporal 
correlations  of  u  are  plotted.  The  abscissa  is  x  for  the  autocorrelation  data 
and  A x/U  for  the  spatial  data.  The  ordinate  is  (1-R).  Cramer's  data  show 
close  agreement  between  spatial  correlations  and  temporal  autocorrelations 
during  both  day  and  night  for  the  60  seconds  of  data  which  were  plotted.  This 
shows  that  Taylor's  hypothesis  is  useful  in  the  atmosphere  as  well  as  in  the 
laboratory. 

Figures  5*20  and  5.21  of  [Lumley  and  Panofsky,  1964]  depict  the 
5Ti  autocorrelation  functions  of  both  u  and  v  lagged  in  both  space  and  time  for  a 

time  period  of  20  seconds.  For  the  daytime  observations  agreement  is  extremely 
zii  close,  but  time  autocorrelations  tend  to  be  slightly  larger  than  space 

\  autocorrelations.  During  a  typical  night  period  space  and  time  autocorrela- 

v;  tions  of  v  are  also  close,  but  time  autocorrelations  of  u  are  consistently 

much  larger  than  space  autocorrelations. 

[Elderkin  and  Powell,  197l]  also  tested  the  space-time  relationships 
in  atmospheric  turbulence.  Their  Table  1  contains  R(x)  and  R(Ax)  for  u' ,  v’ , 
i  and  w’ .  According  to  Taylor's  hypothesis,  R(x)  ■  R(Ax)  if  x  -  Ax/3.  Their 

table  goes  to  252m,  which  corresponds  to  almost  40  seconds  for  the  applicable 
. ,i  3  of  6.4m/sec.  For  u'  the  correlations  R(Ax)  and  R(x)  are  nearly  identical  to 

48m  (7*5  sec),  but  beyond  this  point  some  of  the  R's  differ  oy  20  percent  or 
|  more.  For  v'  and  w'  the  R's  diverge  considerably  after  4  seo.  Beyond  a  few 

1  seconds  all  three  spatial  correlations  are  higher  than  the  corresponding  time 

j  correlations.  Therefore,  one  must  use  some  oaution  when  applying  Taylor's 

widely  used  hypothesis. 
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[Tennekes  and  Lumley,  1972],  in  their  Figure  8. 2,  illustrate 
idealized  behavior  of  R  on  a  spatial  scale*  The  spatial  correlation  of  the  u 
component  deorerses  smoothly  from  unity  at  zero  separation  to  zero  at  large 
separations*  The  autocorrelation  of  the  v  component  decreases  smoothly  to 
eero  and  becomes  slightly  negative  before  leveling  off  at  eero  for  large 
distances* 

Aotual  autocorrelations  may  have  quite  irregular  variations  with 
time.  This  is  illustrated  in  Figures  1  and  2  of  [Stewart,  1975]*  The  figures 
contain  a  representative  sample  of  autocorrelations  which  were  oomputed  for 
1-sec  intervals  from  lag  zero  through  lag  120  for  u,  v,  and  w  at  10m.  In  some 
sets  of  measurements  analyzed  by  Stewart,  the  autocorrelation  function  crossed 
zero  several  times.  On  the  other  hand,  in  some  instances  the  autocorrelations 
of  u  and  v  did  not  reach  eero  during  the  entire  120  sec  for  which  the  com¬ 
putations  were  done.  The  autocorrelation  of  the  w  component  usually  reached 
zero  in  20  seconds  or  less. 

[Mackey  and  Ko,  1975]  show  the  autocorrelation  functions  up  to  150 
sec  for  the  longitudinal  fluctuations  at  13,  28,  43,  and  61m  in  a  typhoon. 
These  curves  decrease  rapidly  for  approximately  25  sec  and  then  fluctuate 
irregularly  about  a  mean  value.  The  61 -m  autocorrelation  levels  off  to  fluc¬ 
tuate  about  a  mean  near  0.3,  and  the  13-m  autocorrelation  levels  off  to  fluc¬ 
tuate  about  0.15*  The  28-m  and  48-m  autocorrelations  fall  between  the  ones 
for  61m  and  13m. 

In  spite  of  the  irregularities  in  many  autocorrelation  functions 
which  are  obtained  from  measurements,  many  investigators  have  attempted  to 
develop  simple  analytical  approximations.  One  of  the  simplest  and  most  widely 
used  approximations  for  the  u  component  is 

R(t)  -  exp  (-./T)  (7) 

where 

T  -  R(T)dx  (8) 

is  an  integral  time  scale*  [Hanna,  1979]  averaged  R(t)  for  the  u  component 
over  several  unstable  runs  made  above  flat  farm  land  in  Minnesota.  This  auto¬ 
correlation  function  was  plotted  versus  %  for  60  sec  and  compared  with  the 
curve  obtained  by  fitting  observations  to  Equation  (7).  This  approximation 
appeared  good  to  within  20  percent  for  the  averaged  unstable  runs. 

[Fichtl  and  McVehil,  1 970]  considered  an  exponential  equation  simi¬ 
lar  to  Equation  (7)  and  applied  it  to  space-lagged  autocorrelations.  They 
examined  a  large  amount  of  data  from  unstable  and  neutral  atmospheres  and 
established  a  dimensionless  length  scale.  For  unstable  conditions  this  scale 
was  within  20  percent,  but  for  neutral  conditions  the  error  was  nearly  45  per¬ 
cent.  Because  high  winds  are  most  frequently  associated  with  a  neutral 
atmosphere,  this  result  would  suggest  caution  in  using  an  exponential  approxi¬ 
mation  for  many  practical  studies. 
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[Maokey  and  Ko,  1975]  fit  their  typhoon  data  whioh  leveled  off 
instead  of  going  to  aero  with  the  more  complicated  expression 

R(*)  ■  a0#“*T  *  «i  oos(mT).  (9) 

They  claimed  that  the  addition  of  the  cosine  term  gave  a  good  representation 
to  their  data. 

[Cramer,  1959]  tried  to  fit  autocorrelations  of  the  u  component  by 
the  equation 


1-R(*)  -  0*2/3  (10) 

where  the  constant  c  is  selected  to  fit  the  data.  Cramer  discovered  that  such 
a  law  fit  some  daytime  experiments  where  the  level  of  turbulence  was  high. 

For  other  case3,  Equation  (10)  we-.s  not  even  approximately  valid  beyond  a  few 
seconds. 

[Frost  et  al.,  1978]  suggest  an  even  more  complicated  function  for 
the  longitudinal  correlation  function 


1>(4l)  ‘  (aLp)^,/,  ^«Lp)  •  (1,) 

where  K  is  the  modified  Bess>3l  function  of  the  second  kind  and  Lp  is  the 
longitudinal  isotropic  turbulence  integral,  scale.  Equation  (ll)  is  referred 
to  as  the  von  Karman  longitudinal  correlation  function,  and  a  simple  exponen¬ 
tial  model  as  the  Dryden  longitudinal  correlation  function.  The  Dryden  func¬ 
tion  is  more  commonly  used  because  there  is  no  compelling  evidence  that  the 
von  Karmen  model  is  better,  and  the  Dryden  function  is  much  simpler. 

IV.  THE  EKMAN  LAYER 

Most  of  the  previous  discussion  has  been  concerned  with  the  lowest 
tens  of  meters  of  the  r in- sphere,  where  the  horizontal  wind  stresses  are 
assumed  to  be  nearly  constant  and  the  wind  does  not  turn  significantly  with 
height.  This  layer  is  sometimes  called  the  surface  boundary  layer,  constant 
stress  laye",  or  constant  flux  layer. 

The  planetary  boundary  layer,  which  is  also  called  the  friction 
layer  or  the  atmospheric  boundary  layer,  extends  from  the  surface  of  the  earth 
to  the  geostrophic  wind  level  [Huschke,  1959].  The  planetary  boundary 
layer  includes  the  surface  boundary  layer  and  the  Flktnan  layer.  Above  the 
geostrophic  wind  level  is  the  free  atmosphere. 

The  Ekman  layer  lies  between  the  surface  boundary  layer  and  the  free 
atmosphere.  An  idealized  mathematical  description  of  the  wind  distribution  in 
this  layer  is  called  the  Ekman  spiral  [Huschke,  1959].  This  Ekman  spiral  is 
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derived  by  assuming  that  within  the  planetary  boundary  layer  the  eddy  viscos¬ 
ity,  K,  and  density,  p,  are  oonstant.  The  motion  is  assumed  to  be  horisontal 
and  steady,  the  isobars  are  straight  and  parallel,  and  the  geostrophio  wind  is 
oonstant  with  height*  The  geostrophio  wind  is  represented  by  the  equation 


where  Ug  is  the  speed  of  the  geostrophio  wind,  f  is  the  Coriolis  parameter, 
p  is  pressure,  and  n  is  horisontal  distance  perpendicular  to  the  flow.  The  n 
axis  increases  to  the  left  of  the  flow  in  the  northern  hemisphere*  If  the  x 
direction  is  now  taken  as  parallel  to  the  isobars  and  positive  in  the  direc¬ 
tion  of  the  geostrophio  wind,  one  can  derive  the  equations  [Hess,  1959] 

u  ■  Ug  (1-e"aB  cos  as)  (13a) 

and  -  - .  .  •  •  . -  - 


v  •  Ug  e”a*  sin  as  (13b) 

where  a  ■  f/2K.  At  a  «  0  the  wind  speed  is  sero.  The  limiting  value  of  the 
angle  of  the  wind  with  the  isobars  as  the  surface  is  approached  from  above  is 
45  degrees,  and  the  wind  points  toward  lower  pressures.  The  wind  vector  turns 
clockwise  with  altitude  in  the  Northern  Hemisphere  and  becomes  parallel  to  the 
isobars  at  the  geostropio  wind  level.  At  this  level,  which  is  near  1  km,  the 
wind  speed  is  slightly  greater  than  the  geostropic  value. 

Another  approach  is  sometimes  used  by  investigators  who  are 
interested  in  levels  above  the  lowest  20  to  40m  if  they  do  not  wish  to  go  much 
above  150m.  [^anofsky,  1973]  shows  that 

u*  -  u*o  -6fs  (14) 

to  a  very  good  approximation.  The  symbol  u#0  represents  the  surface  friction 
velocity.  Under  neutral  conditions  and  homogeneous  terrain  the  following  can 
be  written: 


where  ^  is  the  mean  wind  speed  at  a  reference  level  .  The  exponent  m  has 
been  found  empirically  to  deciease  with  increasing  lapBe  rate.  [Zhang,  198l] 
compared  the  power  law  with  the  simple  logarithmic  wind  law  which  applies  to 
neutral  conditions.  Wind  data  from  a  164-m  tower  were  examined  for  one  year 
in  Nanjing,  China.  For  the  height  range  from  16m  to  164m,  the  power  law 
represented  the  actual  wind  speed  distribution  better  than  the  logarithmic 
law. 


V.  SUMMARY 

Wind  variation  with  height  in  the  surface  boundary  layer  can  be 
approximated  by  a  logarithmic  wind  law.  Thia  law  is  particularly  uBtful  in 
many  applications  because  it  is  quite  good  when  atmospheric  stability  is 
neutral,  and  high  wind  speeds  are  typically  associated  with  neutral 
atmospheric  stability.  When  stability  is  not  neutral,  better  aocuracy  oan  be 
obtained  by  using  an  equation  which  contains  a  small  stability  term  in  addi¬ 
tion  to  the  logarithmic  term.  If  measurements  are  inadequate  to  compute  the 
stability,  one  can  estimate  it  by  Pasquill's  method  which  depends  upon  time  of 
day,  cloud  cover,  and  mean  wind  speed. 

The  intensity  of  turbulence,  ou/0,  varies  in  space  and  time.  It  is 
usually  greater  over  land  than  over  water,  and  the  intensity  is  greater  over 
rough  land  surfaces  than  over  smooth  terrain.  Intensity  of  turbulences  typi¬ 
cally  decreases  rather  rapidly  in  the  lowest  20m  and  decreases  slowly  with 
altitude  above  this  level.  Intensity  of  turbulenoe  normally  is  greater  under 
unstable  conditions  than  under  stable  conditions. 

Some  investigators  prefer  to  measure  intensity  of  turbulence  by  av/C 
instead  of  ou/u.  Near  the  surface  ov/ou  is  less  than  unity  in  slightly 
unstable,  neutral,  and  stable  conditions.  In  moderately  unstable  conditions 
the  ratio  is  near  unity.  As  the  atmosphere  becomes  very  unstable,  ov/ou 
becomes  greater  than  one.  At  higher  altitudes  ov/ou  is  typically  near  unity. 


The  autocorrelation  function  of  u  is  often  irregular  in  individual 
cases,  but  is  somewhat  smoother  when  a  mean  over  a  large  amount  of  data  is 
taken.  A  simple  exponential  function  is  sometimes  used  to  approximate  the 
decay  of  the  autocorrelation  with  time  or  distance,  but  there  is  evidence  that 
this  can  lead  to  errors  of  20  to  45  percent. 

Above  the  surface  boundary  layer  is  the  Ekman  layer  where  the  wind 
approximately  follows  an  Ekman  spiral.  The  planetary  boundary  layer  consists 
of  thia  Ekman  layer  and  the  surface  boundary  layer.  In  the  free  atmosphere 
above  the  planetary  boundary  layer  surface  friction  with  the  earth  has  a 
negligible  influence. 
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